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ABSTRACT 

We performed magnetohydrodynamic simulation of a formation process of 
coronal mass ejections (CMEs), focusing on interaction (reconnection) between 
an ejecting flux rope and its ambient field. We examined three cases with different 
ambient fields: no ambient field, and cases with dipole field of two opposite 
directions which are parallel and anti-parallel to that of the flux rope surface. As 
a result, while the flux rope disappears in the anti-parallel case, in other cases 
the flux ropes can evolve to CMEs and show different amounts of rotation of the 
flux rope. The results imply that the interaction between an ejecting flux rope 
and its ambient field is an important process for determining CME formation 
and CME orientation, and also show that the amount and direction of magnetic 
flux within the flux rope and the ambient field are key parameters for CME 
formation. Especially, the interaction (reconnection) plays a significant role to 
the rotation of the flux rope, with a process similar to "tilting instability" in a 
spheromak-type experiment of laboratory plasma. 
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Introduction 



Coronal mass ejections (CMEs) are one of the most spectacular explosive phenomena 
in which large amounts of mass and magnetic flux are ejected to the interplanetary space. 
CMEs are believed to be consequences of sudden release of magnetic energy, i.e., disruptions 
of coronal magnetic field, as same as flares which are strong brightening observed in soft 
x-ray. The energy released upward can be converted to kinetic energy as ejected flux ropes 
and/or CMEs, while energy released downward can be converted to thermal energy or high 
energy particles. Therefore larger flares sh ow higher association rate with CMEs (IKahler 
I992I : iKim et al.l l2005t lYashiro et al.l 120051 ) , but they do not necessarily accompany each 
other. Association of flares can be explained by energy release rate which depends on the 
condition of coronal plasma. If the energy release is fast enough to heat chromospheric 
plasma, the heated plasma fllls above coronal loops, which are accompanied with strong 
brightening in soft x-ray observed as flares. If the energy release is not fast enough to 
overcome radiative energy loss, the strong brightening does not occur. In the same time, 
if ejecting structures successfully escape to interplanetary, the event can be understood as 
CMEs without flares. On the other hand, in some cases, if an ejection associated with a 
strong soft x-ray brightening, cann ot escape to the inte rplanetary space, the event should be 
observed as a flare without CME (lYashiro et al.ll2006l ). 



The disruptions of coronal magnetic fleld are observed to occur commonly in a wide 
spatial range. Yohkoh satellite reveals that faint giant coronal arcades which are s imilar 
to flares are formed associated with CMEs (IHiei et al.l llQQSt [McAllister et al.lll996[ ). On 



the other hand, the similar evolutions where flares occurs wit h mass (plasmoid) ejections 
are observed accompanied much smaller scale emerging flux ( ISakajiri et al.l 120041 ). Based 
on such observational evidence, IShibatal ( 119991 ) proposed that giant arcades, flares, mass 
ejections, and CMEs are understood with an unifled view so-called 'flares'. The model is 
supported by the statical characteristic of flares. S olar flares show a "power-law" dist r ibution 
betw e en their energ y scales and occurrence rate (jPrakel Il97ll : iDatlowe et al.lll974j : iDennis 
19851 : IShimiz^ll995[ l. 



In contrast to f lares, CMEs show log-normal distribution between the ener gy scales and 
the event number ( lAoki et al.l l2004l : lYurchyshyn et al.l 120051 : iLara et al.l 120061 ) . Statistical 
characteristics of CMEs and flares are different in spite of their common origin. This dif- 
ference in statistical distribution means the existence of some kinds of fllter effects which 
prevent small scale mass ejections to evolve CMEs selectively. Candidates of such fllter 
effects are interactions between eruptions and their ambient magnetic fleld whose spatial 
scale is larger than that of the eruption regions. Important points in formation process of 
CMEs can be not only the formation of ejecting structures (flux rope) but also whether they 
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overcome the obstacles. 



Trigger processes of such eruptions have b e en numerically investigated for a few decade 



(Forbes & PriestI 


1995 


Kliem & Torok 


20061). 



I999I: 


Chen & Shibata 


2000 




Kusano et al. 


2004; 


del ( 


Antiochos et al. 


1999h. the interaction between 



a sheared arcade and its ambient field is considered to play a fundamental role for the trigger- 
ing. In the model, the field configuration is assumed to be favorable to reconnection between 
the erupting sheared arcade and its ambient field. During the eruption the reconnection can 
reduce the poloidal magnetic flux of the flux rope which is supplied by another reconnection 
occurred inside the sheared arcade. As the results of competition between the magnetic flux 
reductio n and the supply , an eruption is capable of f aihng to produce a CME. In some other 



mode ls ( iForbes fc PriestI Il995t iKusano et al.l |2004J : llnoue fc Kusand l2006l : iKliem &: Torok 



20061 ) ideal MHD processes (loss of equilibrium or instability) of detached flux ropes are 
dominant in the trigger processes. In such models, the direction of its ambient field can 
be independent of the trigger processes. If the direction of the ambient field is favorable to 
reconnection, the reconnection can reduce the magnetic flux as mentioned above. In con- 



(Choe & Lee 


1996: 


Shiota et al. 


2008) 



ejecting flux ropes and their ambient field can be a key factor for formation of CMEs. 

From the point of view of space weather science, when CMEs containing southward 
magnetic flux impact the Earth magnetosphere, they make large influences on the space envi- 
ronment near the Earth. Hence, it is also important for space weather forecast to understand 
the whole evolution process of CMEs; how CMEs are formed and how much southward mag- 
netic flux CMEs contain when they reach to the orbit of the Earth. Estimation of southward 
magnetic flux is a complicated task important for space weather, and requires understanding 
of initial conditions and whole evolution processes of eruptions. For example, in nonlinear 
evolution of kink instabi lity of a strongly t wisted flux rope, its u ntwisting (writhe) motion 



can rotate the flux rope (IKliem et al.l 12004 ; iTorok &: Kliemll2005l ). The ejecting twisted flux 



ropes could change their directions due to the writhe rotations since their onset or formation 
if the ropes become unstable in the mode. 



Rotations of ejecting flux ropes are reported in CME observations. lYurchyshynl (|2008[ ) 
investigated the angles of magnetic clouds, halo CME axes derived from coronagraph images, 
and EUV post-eruptive arcade axes in CME source regions. They found that there are good 
correlation of angles between the axis of EUV arcades and CMEs, and that between axes 
of magnetic clouds and coronal neutral lines (Heliospheric current sheets), while there is 
low correlation of angles between CMEs and magnetic clouds. The fact means that some 
fraction of ejected magnetic structures are significantly rotated during the travel from the 
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solar corona to 1 AU. During the way, the ambient field should be of two type of structure; 
closed and open. In the outer corona or interplanetary, all magnetic field should be open due 
to dragging of the solar wind while most of magnetic field lower corona should be closed. The 
rotation ma y resu lt from the interactions with the both type of ambient field. Especially, 



Yurchyshynl (120081 ) reported that a few events show the deviation between the neutral line 
angles of post eruptive arcades (in lower corona) and the axes of CMEs (in a few solar radii). 
This means that the rotation during the evolution in lower corona, i.e., interaction with 
closed field, is not negligible. 

In this study, we performed a three-dimensional MHD simulation of a twisted fiux rope 
ejected from an active region surrounded by a global ambient magnetic field. According 
to the numerical results, we discuss the condition for formation of CMEs and evolutions of 
CME magnetic field structures, focusing on the interactions between ejecting flux ropes and 
their ambient closed field. 

In the next section, we describe the detailed methodology of the numerical simulation. 
In section 3, the numerical results of different ambient field are shown. We discuss the 
quantitative relations suggested by the numerical results in §4. Finally, we summarize this 
paper in §5. 



2. Numerical Model 

2.1. Numerical Scheme 

We performed three dimensional MHD simulation which solves time variation of 8 physi- 
cal quantities; density, gas pressure, velocity and magnetic field {p,p, v, and B, respectively). 
In order to simplify the MHD equations, the 8 physical quantities are normalized with the 
following typical values in the corona: Lq = Rq = 6.99 x 10^ [km], Bq = 30 [Gauss], 
Pq = mH X 10^ = 1.67 X 10-15 [g cm'^], vq = Vo,a = Bo/y/A^ = 2071 [km s-^], as these 
results. To = ta = 350 [s], po = PqVq j^ = 7.16 x 10^ [erg cm"^]. The normalized MHD 
equations (equation of continuity, equation of motion, induction equation, and equation of 
energy with gravity) are expressed as follows: 

^ + V-(pv) = 0, (1) 

do'v 

^ + V ■ (pvv - BB + ptI) = pg, (2) 

— + V ■ (Bv - vB + #) = - V X {r]3) , (3) 
at 
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de 
dt 



+ V - [(e + pr)v- (v-B)B] = 0, 



where 



pv 



J = V X B, 

P B2 Go 

+ + ^ + — P, 
7 — 1 2 r 



Pt = P + 



B^ 



Go 



2 ' 
GMq 



(4) 

(5) 
(6) 

(7) 
(8) 
(9) 



G is gravitational constant, I is the unit matrix in the Cartesian coordinate, if) is additional 
variables described below, and Mq is mass of the Sun. Divergen ce operator is solved by a 
finite volume method with the HLLD nonlinear Riemann solver (IMiyoshi &: Kusand |2005[ ) 
which is combined with the 3rd order Monotone Upstream-centered Schemes for Conserva- 



tion Laws (MUSCL) developed by Ivan Leerl (119791 ) and the 2nd order Runge-Kutta time 



integration. Specific heat ratio 7 is set to be 1.05 which means an additional heat source. 
Resistivity r] is set to be which means reconnection is caused by the numerical diffusion. 

Generally, the results of time integration of equation dH]) multi-dimensional MHD sim- 
ulation is possible to cause numerical V ■ B due to discretization of numerical domain. The 
equations ([2]), (|3]), and (jl]) are derived with Solenoidal condition, and therefore such finite 
numerical V ■ B can break down the simulation. In order to keep the numerical V ■ B to be 
sufficiently small value, we apphed an additional variable ip and additional equation. 



^ + 4V-B + 4/cJ^ = 



(10) 



based on the hyperbolic divergence B cleaning method (jPedner et al.ll2002l ). The coefficient 
Ch in Equation ffTOj) is propagation speed of the numerical V ■ B, and the coefficient Cp i s 
diffusion coefficient of ip. This numerical methods are described more in lShiota et al.l ( 120081 ). 



The numerical domain is a spherical shell of (rmin,^ max) = (1.01,5.0). The domain is 
discretized with Yin- Yang grid ( Kageyama fc Satdl2004j ). a chimera grid composed of two 
congruence spherical partial shells (vr/4 < 9 < 37r/4,— 37r/4 < < 37r/4) each of which is 
installed in different direction. The physical quantities on the edges of the each partial shell 
are substituted by interpolated values at the same position in the counterpart shell. 

Thus the boundary condition in this grid system only at inner and outer edges for the 
radial direction. The inner boundary is assumed to be the line-tying wall, where -Bnorm is 
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constant and v = 0. The outer boundary is assumed to be the free boundary, where the 
radial gradients of all quantities are zero. 

Note that, in this simulation, test particles which are distributed in some manners in 
initial condition are advected with plasma bulk flows. The magnetic fleld lines in Figures in 
this paper are obtained with the integration of the magnetic fleld from the position of the 
test particles of each time step. Thus we capture the evolution of individual magnetic fleld 
lines with this technique. 



2.2. Initial condition 

The subject of the present work is to understand how a flux rope ejected from an active 
region evolves to a CME. The initial magnetic fleld consists of the ejecting flux rope (Bs) and 
ambient global fleld (BD).We performed three cases of simulation where only the ambient 
fleld Bd is different. The ambient global fleld (Bd) is chosen to be simple dipole fleld: 

, ^ , 2Br, sin 6 , , 

i?r,D(r,g,0)= (11) 

ge,D(r,g,0)= . (12) 

In order to study the interaction between an ejected flux rope and ambient fleld, we carried 
out three cases of simulation with different strength and direction of the ambient flelds, 
i.e., their amplitude -Bd of 0.0,— i?D,0! and -Bd,0! which are labeled with case A, B, and C, 
respectively. -Bd,o = 0.124 (3.72 [G]) is a base fleld strength of the corona where the Alfven 
speed is equal to the sound speed of T = 2 x 10^ [K]. 

We mimic the ejected flux rope formed as a result of eruption in an active region with 
a spheromak type magnetic fleld which is a linear force-free fleld in a completely isolated 
sphere (see Figured^). 

Bs(r,^,0) = B(f,^~,0), (13) 

where f = a(r — rs) is the translated and rescaled local spherical coordinate whose origin 
corresponds to the center of the spheromak. The spheromak type fleld is expressed as follows: 

Brij, ^, 0) = - 25so^ COS ~Q, (14) 

S,-(f, ~d, 0) =5so + j( (r)^ sin ^, (15) 

5.(f,^,0)=5soJi(r)sin^~, (16) 
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where Bso is field strength of tlie splieromak and ji andj[ are first order spherical Bessel 
function and its derivative: 

. , , sinx — xcosx 

ji(x)= , (17) 

, 2x cos X — (x^ — 2) sin X , , 

Ji W = -3 • (18) 

We adopt only the field within the isolated spherical shell r < a (transparent sphere in 
Figure [1^) while Bs = outside of the shell. As definition of the spherical Bessel function, 
a = 4.493409458 and a = a/as- The toroidal flux ($s) and poloidal flux ($p) of this field 
are obtained by numerical integrations; 

$s = / d9 fB^df = 0.261alBso, (19) 
Jo Jo 

$P = y d0 y - Bg(9 = -) df = 0.0736alBso. (20) 

In the present paper, we choose the position of the center of the spheromak rs = 
(l.l,7r/2,0) in the spherical coordinate (on the x axis), and the major axis of turns is 
parallel to the z axis. The radius of the spheromak as = 0.2, i.e., 1.4 x 10^° [cm]. Here, we 
assume the toroidal flux of the spheromak $s = 0.0174 = 2.55 x 10^^ [mx] which is a typical 



amount of magnetic flux within a CME (IKataoka et al.ll2009l ). Then Equation flTIJl) yields 



the fleld strength of the spheromak Bso = $s/(0.261ag) = 1.67 = 50 [Gauss]. 

Although a spheromak is a linear force-free fleld, we adopt it only within the boundary 
f < a. At the spherical boundary r = a, the local radial (Bf) and azimuthal (B^) components 
vanish and then only the local zenithal component Bg exists while B = outside of the 
boundary. Therefore the force does not balance at the spherical boundary, i.e., a strong 
outward gradient of magnetic pressure exists there. The Lorentz force due to the current 
which corresponds to a steep magnetic pressure gradient causes the strong expansion. The 
spheromak fleld naturally swells in the simulation. 

Figures [T] shows schematic pictures of magnetic fleld direction on the meridional plane 
through the center of in cases B ([T]d) and C ([It). The spheromak local zenithal fleld on its 
surface r = a satisfles B^ {r = a) > (cf. eq. (fT5|) ). which is southward. On the other hand, 
the ambient fleld is northward {Bd < 0) in case B, while it is southward in cases C. On the 
equator, Bg (f = a) and -Be,D < becomes anti-parallel in case B, while parallel in case C. 
Hence, cases A, B, and C are referred as 'no ambient' 'anti-parallel', and 'parallel' cases in 
this paper. 



Density and pressure are determined by hydrostatic equilibrium with uniform temper- 



- 8 - 



ature T = 2 x 10*^ [K], 



p = Ph exp 
P 



Go7 A _ ^ 
Ph \r 



1- 



p_ 
Ph 
0, 



(21) 

(22) 
(23) 



where pb = '^Po^bT / {mupo) is the normahzed pressure on solar surface (r = 1). In the 
present study, we did not apply the solar wind, which is also important factor for the for- 
mation of a CME. This is because we intend to purely investigate the interaction between 
the ejecting and the ambient fields in the early phase of the CME formation. 



3. Numerical Results 

3.1. Common evolution 

In all cases, the spheromak (flux rope) initially swells outward due to the non-equilibrium 
initial condition and the line-tying boundary as described in the previous section. At the 
same time, the expanding flux rope also shows the similar rotation around the line of the 
ejection in all case. This rotation is caused by an untwisting (writhe) motion of the highly 
twisted initial spheromak structure (see Figure [TK). 

In order to show the height-time evolution, the right panels of Figure [2] illustrate time 
variation of plasma /3 = O.bp/B"^ along the x axis. The low plasma /3 regions along the x 
axis express the positions of the flux rope. Therefore the low plasma (3 in the height-time 
diagrams illustrated in Figure [2] corresponds to the trajectories of the flux ropes. In all the 
cases the flux rope rapidly expands just after the start, followed by rising and expanding 
motion with some deceleration. 

In all the cases, many small features appear around the altitude r ~ 1.1 and rise as time 
below the flux rope in the right panels of Figure O These structures are small plasmoids 
formed by tearing of the current sheet between the anchored axial flux. Once such plasmoids 
are formed in the current sheet, due to gradient of magnetic pressure, they are ejected to 
upper direction, and flnally collide and coalesce into the above flux rope. However, the 
momentum and energy additions by the plasmoids are t oo small to accelera te the flux rope. 



The evolution is similar to the results of 2D simulation ( IShiota et al.ll2008l ). 
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3.2. 'No ambient' case 

In 'no ambient' case (case A), tlie flux rope can expand into unmagnetized plasma. We 
simulated this case as a reference case with no magnetic obstacle. The expanding flux rope 
drags outward its anchored axial flux. The dragged axial flux forms fi-shape and at the 
center of them a current sheet is formed. Although the reconnection occurs between the 
dragged fluxes, a part of the flux remains anchored. 

Figures [2^ and Wp show the magnetic fleld structure at t = 50 (well-developed phase) 
and the height-time diagrams in this case. Figure [2]d expresses the time evolution of the flux 
rope in this case. The flux rope rapidly expands just after the start (t < 10: early phase), 
and then it shows rising and expanding motion during middle phase (10 < 30). The flux 
rope is slightly decelerated. In the late phase {t > 30), the flux rope expands self-similarly 
at almost a constant speed, and flnally reaches the outer boundary of the numerical domain. 



3.3. 'Anti-parallel' case 

We simulated the 'anti-parallel' case (case B) as a condition where the directions of 
poloidal fluxes of the ejection and the ambient fleld are opposite. Such a condition is probable 
if the flux rope is formed and ejected in a newly emerged flux system whose direction is 
opposite to the pre-existing a mbient fleld. This situa tion is modeled as a candidate system 



for rapid eruption for a CME (lAntiochos et al.lll999l ) 



The time evolution of the flux rope is signiflcantly different from that in case A. Because 
the direction of the ambient fleld is opposite to the poloidal field of the ejecting fiux rope, 
th ere is a strong current sheet where magnetic reconnection occurs. Similar to the scenario 



of lAntiochos et al.l ( 119991 ) . the rise motion of the flux rope presses the current sheet, and 
then enhances its current density, i.e., the reconnection rate there. However, the poloidal 
flux of the flux rope is peeling off due to the reconnection with the ambient as it rises. 
Finally, as shown in Figure after most of the poloidal flux is reconnected, the flux rope 
no longer exists. The magnetic flux becomes just coronal loops within which shear Alfven 
waves propagate. Such a result may be a candidate condition for plasmoid ejections (flares) 
without CMEs. 

We can see the time evolution of the flux rope in height-time diagram (Figure |2jl). In 
the early phase {t < 10), the trajectory and the vertical size are similar to those in case A. 
After the phase, the top edge of the flux rope becomes lower than that in case A because 
of the flux peeling (10 ^t<35: middle phase). Finally (t > 35: late phase), the plasma (3 
inside the flux rope becomes comparable to unity, it stays almost the same height. In this 
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phase, the structure is recognized as twisted loops rather than a flux rope. 



3.4. 'Parallel' case 

We simulated the 'parallel' case (case C) as a condition where the directions of poloidal 
fluxes of the ejection and the ambient field a re parallel. CME fo rmation process in a sim- 



ilar but axis ymmetric sit u ation is studied by IShiota et al.l (|2008[ ). The present study is an 



extension of IShiota et al.l (|2008[ ) to study three-dimensionality of CME formation process. 



The time evolution of the flux rope is a little different from that in case A, as shown in 
Figure [2^ and[2]F. In the early phase (t < 10), the flux rope rapidly expands, and shows rising 
and expanding motion during in the middle phase (10 ^ t < 45), and then a self-similar 
evolution in the late phase {t > 35). The expanding flux rope forms an f2-shape loop and its 
central current sheet where magnetic reconnection occurs, and finally the flux rope reaches 
the outer boundary, similar to in case A. In contrast to the result of case A, the ejecting 
flux rope is significantly deformed (Figures and [2^), although the flux rope is successfully 
ejected to the outer boundary in both cases. The difference should be caused by reconnection 
between the flux rope and the ambient magnetic field. The height evolution of the flux rope 
is similar to that in case A but the flux rope continues to be decelerated (Figures |2]b and[2f). 

We investigated the time evolution of the magnetic field configuration. Figure |3] and H] 
shows the time evolution of the magnetic field configuration seen from the side and the face 
of the ejecting flux rope. Some characteristic field lines are highlighted with colors (green, 
light blue, blue, and red) in the figures. These lines are drawn by integrating magnetic field 
from the same points on the inner boundary in each time step. 

The "green" field line represents the initial axial field line. As the flux rope expands, 
the axial field lines are dragged and stretched as shown in Figure |3^ - 13^. Finally the field 
line reconnects at the f2-shaped center as described in §3.21 

The "light blue" field line indicates "interchanged field line". The line is initially an 
ambient field line (Figure [3^ and UK). As the flux rope expands, the line is pushed aside 
by the flux rope (Figure [3]d andHb). Because the direction of magnetic field is parallel, the 
ambient field line does not reconnect with the flux rope. However, as the flux rope evolves, it 
rotates due to writhe motion. The field lines are connected to inside of the flux rope (Figure 
|3t andHt). The result implies that the field lines reconnect with the flux rope at the current 
sheets formed on the flux rope surface. The location of the reconnection is examined later 

(M- 
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The "blue" field line indicates also "interchanged field line" . The line is also initially an 
ambient field line further than "light blue" line. The "blue" field line reconnects with the 
ejecting fiux rope at the current sheet formed on its surface as same as the "light blue" field 
line. 

Reconnection occurs also in the front of the fiux rope. In the initial condition, the 
poloidal field of the surface of the fiux rope is parallel to the ambient field. As mentioned in 
§3.H however, the fiux rope rotates due to its writhe motion. The direction of magnetic field 
on the fiux rope front rotates counterclockwise as shown in Figure IH The rotation cause the 
formation of a current sheet on the front and then magnetic reconnection occurs there. 

See the "red" field line in Figure [3] and |H The "red" line was initially ambient dipole 
field but finally connected to the inside of the fiux rope, i.e., becomes the new axial field of 
the ejecting fiux rope (Figure [3]F andSf). 

4. Discussion 

4.1. Magnetic Flux relation for CME formation 

The amounts of magnetic fiux and magnetic helicity are key parameters to understand 
the numerical results. Here we estimate the amounts of magnetic fiux and helicity and discuss 
the relation between them and a condition for CME formation. In the initial conditions of 
case B and C, though the spheromak field structure is superposed into the ambient dipole 
field, the ambient dipole field component is so weak that it does not change the spheromak 
field. Hence we can regard the time evolutions in case B and C as the results from the 
interaction (reconnection) between ejecting fiux rope and the ambient field. 

Whether an ejecting fiux rope can evolve to a CME depends on the magnetic fiux 
relation of the overlying ambient field and the ejecting fiux rope. In the present study, the 
dipole field is assumed as the ambient field and therefore we can integrate the amount of the 
fiux with the following equation: 




Because only the magnetic fiux above the spheromak field is relevant to the CME condition, 
the azimuthal width 0min, 0max of the overlying magnetic fiux are determined by that of 
the spheromak field in the initial condition: 0min, 0max = ± tan~^(as/|rs|). Then we obtain 
the overlying magnetic fiux $d = 1.19 x 10^^. On the other hand we can also integrate 
the poloidal fiux of the spheromak field with the equation fl20|) . The azimuthal width in 
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spheromak coordinate is set to ±7r/2 when the upper half of the spheromak field is assumed 
to commit to their ejection. Then we get = 2.45 x 10^'^. The relation between amount of 
magnetic flux $p < $0 in case B and C. This condition naturally explains the disappearing 
of the flux rope due to flux peeling-off reconnection in case B. In case C, the reconnection is 
not completely anti-parallel and so inefficient that the amount of the stripped magnetic flux 
is less than that of the flux rope. In case A, there is no obstacle ambient field. Hence, the 
flux rope in case A and C is not dissipated. 

Although only three typical cases are modeled in the present study, there can be any 
intermediate situations. For example, the ambient field is anti-parallel but total amount 
of magnetic flux is smaller than that of an ejecting flux rope, the flux rope also retains 
its shape and escapes to the outer space. Such a condition is a ca ndidate conditiori wher e 
explosive energy release h as been studied as the break out model (lAntiochos et al.l Il999l ) . 
In the case of small scale, Galsgaard et al. ( 2007 ) and Archontis fc Tor5^ ( 20081 ) studied 
emerging flux into preexisting and weak coronal field. iGalsgaard et al.l (120071 ) showed rising 
motion of the emergin g flux is not affected by th e orientation between the emerging flux and 
the preexisting field. lArchontis fc Torokl (120081 ) showed in anti-parallel case the emerging 
flux rope successfully erupts while remains confined in other case. In such cases of emerging 
flux, the amount of emerged magnetic flux becomes more than that of the coronal field 
within the same spatial scale, because the strength of the magnetic field is very strong in 
the lower atmosphere. However, since the magnetic field in the real Sun is very complicated, 
magnetic field in small scale does not necessarily have same direction with global scale field. 
Furthermore the spatial scale of emerging flux should be (sometimes very) small compared 
to the solar radius, the flux of global scale ambient field is possibly comparable or larger 
than that of flux ropes. 

In this study, we assume a spheromak field structure smaller than the global scale (the 
solar radius) and which has a typical amount of magnetic flux 2.553 x 10^^ [mx]. The ambient 
field ~ 2.5 [G] is also comparable to or a little weaker than a typical coronal value. The 
difference in the direction and the amount of the overlying field leads to much different 
results. This result may give a hint to understar id the observationa l result that all X-class 
flares are not necessarily associated with CMEs (lYashiro et al.ll2005l ). 



Hence, the overlying field above an eruption region is an important factor for CME 
formation, and therefore its detail estimation with a realistic corona model is required for 
CME prediction for space weather forecast. Note that we also did not consider the effect of 
solar wind which makes substantial amount of magnetic flux open and decreases magnetic 
obstacle. The detail estimation of amount of magnetic flux using observational data and a 
realistic corona model is an important task for a CME prediction. 



- 13 - 



Shiota et al.l ( 120081 ) discussed a condition for CME formation with 2D simulation of 
continuously sheared arcade which results in a CME or a confined flux rope formation. They 
suggested that the condition for CME formation depends on the balance between amounts 
of the magnetic helicity (H) within the flux rope and sheared arcade, and overlying magnetic 
flux ($ov) which confines the sheared structure, as follows; 

wi < ^^^^ 

In these symmetric condition, maximum energy for the flux rope to escape requires very 
much due to the work to drag the overlying flux. Therefore the condition of equation fl2S]) 
is possibly a upper limit for a 3D condition. 

The relative magnetic helicity of the each case of the present study, is estimated from 
the initial condition 

= y (A + Ap) ■ (B - Bp) dV^, (26) 

where Bp and Ap are the potential magnetic fleld and the corresponding vector potential 
which are determined by the magnetic fleld radial component on the inner boundary of the 
numerical domain. The estimated relative helicity is H = —3.65 x 10~^, —2.96 x 10~^, —4.34 x 
10~^ in case A, B, and C, respectively. The relation of equation (1251) for case becomes 

\H\ \H\ 4.34 X 10-4 ^ ' 

This relation suggests that the conflnement due to the overlying fleld is inefficient to preven t 



a flux rope escape in case C the results are consistent with the results of IShiota et al.l ( 120081 ). 
As same as discussed above, the solar wind decreases the amount of the overlying magnetic 
flux and help the flux rope to escape from the Sun. 



4.2. Rotation of ejecting flux rope 

As described in ^ the rotations of ejecting flux ropes in the early phase in any cases are 
caused by the writhe motion of the spheromak structure in the initial conditions. However, 
the flnal angles of the escaping flux ropes are quite different (~ 90°) between in case A and 
C (Figure [2^ and [2^). The difference is caused by only the existence of the ambient fleld, i.e., 
the interaction between the ejecting flux rope and the ambient fleld as mentioned in §3.41 
The interaction is the reconnection of twisted fleld within the flux rope to the untwisted 
ambient fleld. As a result of reconnection of not completely anti-parallel fluxes, strongly 
bent magnetic flux can be formed. As the result of a relaxation process of the flux, the flux 
rope can be rotated. 
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Angular momentum [L^) and torque(r^) around the x axis are obtained from their 
integral in the space Vi = {r| rmin r < Tmax, tt/A < 9 < 37r/4, — 7r/4 < (p < 7r/4} 

dy (r X pv) ■ X (28) 





dy (r X F) ■ X (29) 

where F expresses any force of Lorentz force, or gradients of gas pressure, or dynamic pres- 
sure. The time variations of the torques around the x axis are shown in Figure [5^ and that 
of angular momentum x component is shown in Figure [5)d. Figures [5^ and [5)d show that 
the structure continues to be rotated dominantly by the torque due to the Lorentz force. 
This result appears that the rotation caused by the writhe motion of the initial flux rope. 
However, Figure shows the three dimensional distribution of the torque. The net torque 
mainly works at the surface of the flux rope (indicated pink surfaces emphasized with a 
white circle in Figure (St). It appears that the rotation are caused by interaction between 
magnetic field inside the flux rope and the ambient magnetic field. Magnetic field line which 
penetrates the strong torque region is drawn with "purple" line in Figure Et. The "purple" 
line appears strongly wind, i.e., just after reconnected. This result means that the strong 
torque regions correspond to the current sheets where the reconnection between the flux rope 
and the ambient field occurs. 

Figure Oi shows a schematic picture of the evolution of the magnetic field lines in case 
C mentioned above and in §3.41 The expanding flux rope initially rotates counter-clockwise 
due to its writhe motion. As a result, current sheets are formed at two point on the surface, 
where the anti-parallel field are contact (indicated with red lines in Figure [5li). Reconnection 
occurs in the current sheets and then strongly wind field lines (indicated with green and light 
blue lines in Figure 01) are formed. The untwisting motion of the winding field causes a 
strong torque (purple arrows in Figure [Sji) which rotates the flux rope. 

In a few decades ago, a magnetic field configuration of spheromak in equilibrium state 
under uniform field has been applied to a laboratory plasma experiment, where plasma 
must be confined within a device. H owever, such a sphero mak type field configuration is 



unstable to tilting mode perturbation (ISato fc Hayashilll983l ). The field configuration settles 



down where magnetic moment of the spheromak is parallel to the background field. In the 
nonlinear stage of the instability, the magnetic field causes reconnection which is the same 
as explained in Figure EJi. The numerical results in this study imply that such an effect in 
laboratory plasma works also in the coronal plasma. 



Such rotations of ejecting flux ropes are reported in some CME observations. lYurchyshyn 



(120081 ) investigated the axis angles of magnetic clouds, halo CME axes derived from coro- 
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nagraph images, and EUV post-eruptive arcade axes in CME source regions. They found 
that there are good correlation between the axes of EUV arcades and CME angles within a 
deviation of ±45°. The fact means that some fraction of ejected magnetic structures are sig- 
nificantly rotated during the travels in the lower solar corona but the amount of the rotation 
is not so much. 

The numerical results of the present work are consistent with the observation in that 
flux ropes can be rotated significantly by the interaction with closed f i eld. However, the 
rotation angle in case C is more than that of most CMEs in lYurchyshyru (120081 ). The reason 
is that the effect of the ambient field in case C is much stronger than the real solar corona 
because of the following reasons: The dipole field assumed in this study decreases in with 
the distance from the Sun r, and is the potential field whose influence can reach most far 
away. Most of the dipole field flux are closed because solar wind is not taken into account. 
Alfven speed in the simulation is very small because of the assumption of the high plasma 
density. The slow evolutions allow that the ejecting structure can show enough rotation. 
The condition in case C is an ideal case to examine the interaction between the ejecting flux 
rope and the ambient field. Furthermore, although the time evolution in case C is very slow 
(the terminal CME speed ~ 40km s~^), the terminal speed is comparable to Alfven speed 
inside the flux rope. 



The r otation proce s s of ej e cting flux ropes is als o numerically studied by lGibson &: Fan 



(2008 


) and 


Lvnch et al. 


(2009 


)• 



numerical result of eruption of a twisted flux rope within a potential arcade. They showed 
that t he rotation of the flux rope axis (~ 115°) is caused by the writhe motion. iLynch et al. 



( 12 0091 ) discussed the rotation precess with the results of two eruptions of arcades sheared 
oppositely each other. They showed that the opposite chirality of the sheared arcades leads 
the same amount of opposite rotations of the resulting flux ropes. In contrast to these 
studies, the results of the present study show that the rotation of the flux rope is affected 
by not only the writhe motion but the interaction with ambient field as described in ^ 
Contrast to their works, the results of the present study suggest that an ejecting flux rope 
can be significantly rotated as a result of reconnection with its ambient field. 



4.3. Evolution of Anchored Axial flux 

In case C, because of the interchange reconnection mentioned in §3.41 the foot points 
of the ejecting flux rope appear to move outside (colored lines in Figure |3] and S]). In 
the initial condition only the "green" line was included in the anchored axial flux of the 
ejecting flux rope, but finally all the "light blue" , "blue" and "red" lines are included. These 
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changes of anchored foot points are caused by a topological change of magnetic flux due to 
reconnection in two regions; at the surface and behind of the ejecting flux rope. As described 
in §31 the ejecting and expanding flux rope forms fi-type loops of the initial axial flux. The 
first reconnection occurs in the current sheet at the center of the waist of the loop (see the 
"green" line in Figure |3] and H]) . This reconnection severs the anchored axial flux and makes 
the ejecting flux rope structure detached. 

At the same time, the ejecting flux rope reconnects with the ambient field in its side 
("light blue" and "blue" lines) and front (the "red" line). In the initial condition of this 
case, the poloidal field of the flux rope is parallel and therefore unfavorable to reconnection. 
However the flux rope rotates due to its untwisting (writhe) motion, as described in §4.21 
The reconnected ambient field is connected to the inside of the flux rope, i.e., becomes the 
new axial flux of the ejecting flux rope. The interchange of the anchored axial field makes 
distance between the legs larger and makes the foot point field strength weaker. As a result, 
the magnetic tension force can decrease and then the flux rope becomes easier to escape 
from the Sun. 

The interchange of anchored flux explains the evolution of dimming region observed in 
many flares and/or CMEs with soft X-ray or EUV. In many flares, strong dimming are often 
observed in both sides of the post eruptive arca des, and there are following weak dimming 
far from the flare cite (ISterling fc Hudson! 119971 ). Such strong dimming can be formed by 
a sudden expansion of the erupting core and upward dragg ed overlying flux coupled with 
coincident rarefaction due to reconnection (jShiota et al.ll2005l ). The following weak dimming 
are formed by some interaction (stretched or reconnected) with the er upting structure. Thi s 
evolution of the magnetic field structure is similar to that mentioned by lGibson fc Fan! (120081 ). 
The results of case C in the present study also show that the weak dimming are seen in new 
footpoints of axial flux formed due to reconnection described in this section. 



5. Summary 

In this study, we performed three global MHD simulations of interaction (magnetic 
reconnection) between an ejecting flux rope and its ambient field, as a launch process of 
coronal mass ejections (CMEs). The initial conditions are assumed as ejecting flux ropes 
under the three different ambient fields: A case without ambient field and cases with dipole 
field of two opposite directions which are parallel and anti-parallel to that of the flux rope 
surface. The flux rope can escape to the outer boundary of the numerical domain (evolve to 
CMEs) in the no ambient case and one dipole case where the direction of ambient field is 
parallel, while it cannot escape in the other dipole case. 
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The flux ropes show rotations perpendicular to the ejecting direction, that are caused 
by not only untwisting (writhe) motion of the flux rope but also the interaction between 
the ejecting flux rope and the ambient field. Especially, the results show the possibility 
of rotation due to the relaxation of a complicated field structure which results from the 
reconnection between the flux rope and the ambient field. The reconnection and resulting 
rotation are similar to that seen in a spheromak experiment of laboratory plasma. 

The results of this study show that the interaction between an ejecting flux rope and its 
ambient field is a significant process for determining CME formation and CME orientation. 
For the achievement of space weather forecast it is necessary to predict occurrence of a CME 
when an flare occurs. A powerful tool for such prediction is MHD simulation of the flare and 
accompanying ejecting flux rope under a realistic coronal magnetic field. 
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(t)) case B 



(c) case C 
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Fig. 1. — (a) Three dimensional view of initial magnetic field in case C. The surface color 
shows Br map and tubes are magnetic field lines. Transparent surface illustrates the surface 
boundary of the spheromak field, (b) Schematic picture of initial magnetic configuration in 
the y = plane in case B. Spheromak type field of whose radius is as are assumed in the 
initial condition (shadowed region). The solid lines express magnetic field lines projected 
on the y = plane, (c) Schematic picture of initial magnetic configuration in case C shown 
with the same way as panel (b). 
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Fig. 2. — (a) Field configurations in late stages {t = 50) in case A (no ambient), (b) timeslice 
of plasma /3 on the x axis in case A (no ambient), (c) and (d) are those in case B (anti- 
parallel), and (e) and (f) are those in case C (parallel), respectively. In the timeslice images, 
the trajectories of the flux rope are recognized as low /3 region with rising motion and white 
dashed lines indicate timing when the evolution of the flux rope changes describes in §3. 
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Fig. 3. — Time evolution of magnetic field and velocity in case C. Panels (a) - (f) show those 
at t =0, 2, 5, 10, 20, and 40, respectively. The inner boundary is shown with the sphere in 
the left part of each panel, whose colors indicate Br map. The same field lines are displayed 
with tubes in each panel. Characteristic field lines are colored with green, light blue, blue, 
and red. The velocity map at y = are displayed with colors on the vertical transparent 
background. 
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Fig. 4. — Time evolution of face-on view of magnetic field structure in case C. Panels (a) - 
(f) show those at t =0, 2, 5, 10, 20, and 40, respectively. The meaning of the colors on the 
sphere and the tubes are the same as those in Figure [3] 
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■ a) Torque aruund k-^kIs 



(c) 3D tJistribuHon of Torque anouncJ if-axis 




Fig. 5. — Time variation of torques (a) and angular momentum (b) around x-axis in case 
C (parallel). Green, blue, red, and black lines in panel (a) show time variations of torques 
caused by dynamic pressure, gas pressure gradient, Lorentz force, and total force, respec- 
tively, (c) Three dimensional distribution of the magnetic torque around the x axis. Pink 
surfaces show positions where counter-clockwise (negative x) torques work. White circle 
indicates the position of the net counter-clockwise torque which substantially rotates the 
flux rope, (d) Schematic pictures of magnetic field lines on a horizontal cross section in the 
ejecting structure. Due to its writhe motion, the twisted flux rope rotates counter-clockwise 
(evolution form the left panel to the center panel). As a result, current sheets (indicated with 
red lines in the center panel) are formed on the partial surface of the flux rope. Reconnection 
in the sheets produce steep bend magnetic field (green and light blue lines in the right panel 
of (d) and purple tube in panel (c)), which causes the torque due to the magnetic tension 
force. 



